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O ■ 1 Problem 
o 

^ ' Gaussian beams provide the simplest mathematical description of the essential features of a 

Qj • focused optical beam, by ignoring higher-order effects induced by apertures elsewhere in the 

Q . system. 

: Wavefunctions ^(x, t) = ^(x)e-''^* for Gaussian laser beams 01211211110011111011111112] 

of angular frequency u are typically deduced in the paraxial approximation, meaning that 
^ I in the far zone the functions are accurate only for angles 6 with respect to the beam axis 

that are at most a few times the characteristic diffraction angle 



J: = ^ = = 

O ' where A is the wavelength, /c = ci;/c=27r/Ais the wave number, c is the speed of light, Wq 

t^. is the radius of the beam waist, and zq is the depth of focus, also called the Rayleigh range, 

r-| '■ which is related by 

^ . Since the angle with respect to the beam axis has unique meaning only up to a value of 7r/2, 

^ I the paraxial approximation implies that 6*0 -C 1, and consequently that Zq ^ Wq A. 

I The question arises whether there are any "exact" solutions to the free-space wave equa- 

CN ■ tion 

I for which the paraxial wavefunctions are a suitable approximation. For monochromatic 
waves, it suffices to seek "exact" solutions to the Helmholtz wave equation, 

VV + = 0. (4) 

This equation is known to be separable in 11 coordinate systems ^31 which oblate 

^ ! spheroidal coordinates are well matched to the geometry of laser beams, as shown in Fig. ^ 

^ I "Exact" solutions to the Helmholtz equation in oblate spheroidal coordinates were devel- 

oped in the 1930's, and are summarized in ^3^1 EI- These solutions are, however, rather 
intricate and were almost forgotten at the time of the invention of the laser in 1960 ^Hj- 

This problem does not explore the "exact" solutions, but rather asks you to develop 
a systematic set of approximate solutions to the Helmholtz equation in oblate spheroidal 
coordinates, which will turn out to be one representation of paraxial Gaussian laser beams. 
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Figure 1: The x-z plane of an oblate spheroidal coordinate system (C, ^,0) 
based on hyperboloids and ellipsoids of revolution about the z axis, with foci 
at {x,z) = (±2o,0). The coordinates have ranges < C < oo, — 1 < ^ < 1, 
and < < 27r). 



The relation between rectangular coordinates (x, y, z) and oblate spheroidal coordinates^ 

(c,e,</') is 

X = Zq\Ii + C^y^l - ^2 COS0, (5) 

y = ^o\/l + CVl-e'sin0, (6) 
z = ZoC^, (7) 

where the length zq is the distance from the origin to one of the foci of the ellipses and 
hyperbolae whose surfaces of revolution about the z axis are surfaces of constant C and 
^. Coordinate is the usual azimuthal angle measured in the x-y plane. For large C, the 
oblate spheroidal coordinates are essentially identical to spherical coordinates (r, 6, 0) with 
the identification C = r/zo and = cos 6*. 

An obvious consequence of the definitions ©-(IZ|) is that 

r± = ^x^ + = zo\/l + CVl-e- (8) 

It is clear that the oblate spheroidal wave functions will have the mathematical restriction 
that the entire wave crosses the plane 2 = within an iris of radius Zq, the length used in the 
definitions ®-(IZI) of the oblate spheroidal coordinates. In effect, the plane 2; = is perfectly 
absorbing except for the iris of radius Zq. 

You will find that the length zq also has the physical significance of the Rayleigh range, 
which concept is usually associated with longitudinal rather than tranverse behavior of the 

^Oblate spheroidal coordinates are sometimes written with C = sinhu and ^ = cosv or sinw. 
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waves. Since the paraxial approximation that you will explore is valid only when the beam 
waist Wq is small compared to the Rayleigh range, i.e., when wq <^ zq, the paraxial wave 
functions are not accurate descriptions of waves of extremely short focal length, even though 
they will be formally defined for any value of wo. 

The wave equation (jH) is separable in oblate spheroidal coordinates, with the form 

+ ^ )lv + 1^(1 - OlJF + . ^2x., 7:^1^ + k Zq{Q + { )^/; = 0. (9) 



dC ' dC dC 'di (i + c2)(i_^ 

It is helpful to express the wave functions in radial and transverse coordinates that are 
scaled by the Rayleigh range zq and by the diffraction angle ^o, respectively. The oblate 
spheroidal coordinate C already has this desirable property for large values. However, the 
coordinate ^ is usefully replaced by 

-=i# = 4(i-a = ^(i-a = ^Tf^, (10) 

^0 ^0 ^ Wq{1 + C ) 

which obeys a ^ {9/9qY for large r and small 9, and a ~ (r_|_/u;o)^ near the beam waist 
where C 0. 

To replace ^ by cr in the Helmholtz equation Q, note that 2^(i^ = —9lda. In the paraxial 
approximation, ^ ~ 1 (which implies that your solution will be restricted to waves in the 
hemisphere z >0), you may to suppose that 

^ --^da. (11) 

Find an orthogonal set of waves, 

= Z^{C)SZia)e^^^\ (12) 

which satisfy the Helmholtz equation in the paraxial approximation. You may anticipate 
that the "angular" functions Sl^{a) are modulated Gaussians, containing a factor a^^'^e"" . 
The "radial" functions are modulated spherical waves in the far zone, with a leading 
factor e*'^^, and it suffices to keep terms in the remaining factor that are lowest order in the 
small quantity 9q. 

Vector electromagnetic waves E = E(x)e~*'^* and B = B(x)e~*'^* that satisfy Maxwell's 
equations in free space can be generated from the scalar wave functions ^/'^ by supposing the 
vector potential A has Cartesian components (for which (V^A)j = V^Aj ^3]) given by one 
or more of the scalar waves -j/^^e"*^*. For these waves, the fourth Maxwell equation in free 
space is cV x B = dE/dt = —ikE (Gaussian units), so both fields E and B can de derived 
from the vector potential A according to, 

E = |-V X B = zA;A + |-V(V- A), B = V x A, (13) 

since the vector potential obeys the Helmholtz equation Q. 

Calculate the ratio of the angular momentum density of the wave in the far zone to its 
energy density to show that quanta of these waves (photons with intrinsic spin S = 1) carry 
orbital angular momentum in addition to the intrinsic spin. Show also that lines of the 
Poynting fiux form spirals on a cone in the far zone. 
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2 Solution 



2.1 The Paraxial Gaussian-Laguerre Wave Functions 

Using the approximation (fTT|) when replacing variable ^ by a, the Helmholtz equation Q 
becomes 

d .2,diP 4 d dtp l + C^-e^ad^^ 4,^ .2 ^2 N , n 

ac^^ + ^ ) + ^^^^ + (i + c>g. ^ ^ + ^ - ^ °- ^''^ 

This equation admits separated solutions of the form (fT^ for any integer m. Inserting this 
in eq. ()14p and dividing by ip, we find 

19,^ ^2,aZ 4 9 95 + 4,^ .2, 4a ^ 

The functions Z and S will be the same for integers m and — m, so henceforth we consider 
m to be non- negative, and write the azimuthal functions as e^*"^*^. With Am as the second 
separation constant, the C and a differential equations are 

d , ^2^dZ / 4 , ^2n ""^^ \ 
ddS el( Aa 

The hint is that the wave functions have Gaussian transverse dependence, which implies 
that the "angular" function S{a) contains a factor e~°^ = e~^^^^oi'^+'' ). We therefore write 
5" = e~°"T, and eq. (|T7jl becomes 

d^T , ,dT ( ^2 ^ ^ 

The function T(o") cannot be represented as a polynomial, but (like the radial Shrodinger 
equation) this can be accomplished after a factor of extracted. That is, we write 

T = a'^l'^L, oi S = (T™/2g-(T^^ gQ ^^^^ (jj^j becomes 

(i^L , ^dL , , 

a-— + {m + l-2a)— + uL = 0, 19 
acr^ da 

where 

^^=jAm-m-l. (20) 

If z/ = 2n for integer n > 0, this is the differential equation for generalized Laguerre polyno- 
mials L^{2a) [in], where 

, , , -A ( — l)''x'' nx n(n — l)x'^ , . 

^ ' '^^{m + k)\{n-k)\k\ m + 1 2(m+l)(m + 2) ^ ' 
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By direct calculation from eq. (jl9p with v = 2n, we readily verify that the low-order solutions 
are 

L;r = l, LT(2a) = l- L^(2a) = l-^^ + - ^ -. (22) 

° ' ^ m+1' 2^ ^ m+1 (m + l)(m + 2) ^ ^ 

The Laguerre polynomials are normalized to 1 at x = 0, and obey the orthogonality 
relation 

L^ix)L^>{x)x'^e-^dz = ; • -. Snn'. (23) 
[m + ny. 

The "angular" functions 5"™ (a) are thus given by 

5„"^(a)=a-/2e-L-(2cr), (24) 
which obey the orthogonality relation 



Jo n \ I n \ I 2'"+i 7o (m + n)!2™+i ^ ^ 

In the present application, < o" < 1/6*0, on which interval the functions 5*™ are only ap- 
proximately orthogonal. Because of the exponential damping of the 5™, their orthogonality 
is nearly exact for 6*0 ^ 1/2. 

We now turn to the "radial" functions ^™(C) which obey the differential equation |TB|) 
with separation constant given by 



= 4(2ri + m + l), (26) 

On 



using eq. (j^UI) with v = 2n. For large r the radial functions are essentially spherical waves, 
and hence have leading dependence e*'^'^. For small polar angles, where ^ ~ 1, the relation 
implies that r ^ zqC, and kr ^ kz^C = 2C/91, recalling eq. Hence we expect the 
radial functions to have the form^ 

Z(C) =e2*'/^oF(C). (27) 
Inserting this in eq. (jl6p . we find that function F obeys the second-order differential equation 
^2\(d'^F AidF\ JdF 2iF\ /4, , \ 

In the paraxial approximation, 9q is small, so we keep only those terms in eq. ()28|) that 
vary as I/^Q' which yields the first-order differential equation, 

dF 

l + C) ^ = -(C + z(2n + m + l))F. (29) 



^It turns out not to be useful to extract a factor e^^^ jr from the radial functions, although these functions 
will have this form asymptotically. 
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For m = n = we write Fq = /, in which case eq. (|29|) reduces to 

(i + c')| = (C + 0(C-o| = -(C + ^)/, (30) 

or 

'4 = -^. (31) 

This integrates to In / = InC — In (C — i). We define /(O) = 1, so that C = —i and 

/ = = = , ■ (32) 

At large C, / ~ 1/C oc 1/r, as expected in the far zone for waves that have a narrow waist 
at 2; = 0. Indeed, we expect that oc 1/C at large C for all m and n. This suggests that 
differs from / by only a phase change. A suitable form is 



f: 



(33) 

Inserting this hypothesis in the differential equation (j2H), we find that it is satisfied provided 

flm.n = 2?2 + m + 1. (34) 

Thus, the radial function is 

i[kzo'^-{2n+m+l) tan^^ C] 

Z™(C) = e'^^^^^Fr = ^= , (35) 

and the paraxial Gaussian-Laguerre wave functions are 

^m./2 „— (7 r r?i/o^'\pi[fc2o''~'^t— (2n+m.+l) tan"'- ?±m<^] 

V;™(a, 0, C, t) = Z™^™e±^'"'^e-^'^* = ' ^ = . (36) 

The factor e-«(2'T-+"i+i)tan u j^,^ ^j^^ wave functions implies a phase shift of (2n+m + l)7r/2 
between the focal plane and the far field, as first noticed by Guoy [201 for whom this effect is 
named. Even the lowest mode, with m = n = 0, has a Guoy phase shift of 7r/2. This phase 
shift is an essential difference between a plane wave and a wave that is asymptotically plane 
but which has emerged from a focal region. The existence of this phase shift can be deduced 
from an elementary argument that applies Faraday's law to wave propagation through an 
aperture [21], as well as by arguments based on the Kirchhoff diffraction integral ^U] as were 
used by Guoy. 

It is useful to relate the coordinates a and C to those of a cylindrical coordinate system 
{r±,(f),z), in the paraxial approximation that ^ ~ 1. For this, we recall from eqs. ((Tj), ^ 
and (Uni) that 

^ = l-9la^ 1, (37) 



6 



so 



and hence 



^ =7271 ^72(1 + ^oM, (3^ 



v2\ / ^2' 



1 + C) ^ wgor 1 + ^ + ^ w^a 1 + ^ ) , (39) 



\ ^0 / \ ^0 / 

where we neglect terms in in the lowest-order paraxial approximation. Then, 

C^±fi + gg^U£fi+ « --l ^ (41) 



For large z eq. (pT|l becomes 

















2^2 1 
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(42) 



as expected. That is, the factor q^{^^(>^-'^^) in the wave functions ()3fjj) implies that they are 
spherical waves in the far zone. 

The characteristic transverse extent of the waves at position z is sometimes called w{^z\ 
From eq. (j40|) we see that the Gaussian behavior e~" of the angular functions implies that 



w{z)=wJ\^''—. (43) 
V ^0 

The paraxial approximation is often taken to mean that variable C is simply zj z^) every- 
where in eq. fj36p except in the phase factor e*'^^°'', where the form (PT|) is required so that 
the waves are spherical in the far zone. In this convention, we can write 



m/2 -o- T m Ir, \ {{kzXl+r"^ /2(z2+22)]-wt-(2n+m.+l) tan-i(z/2o)±m</>} 

0, t) = ^— = . (44) 



The wave functions may be written in a slightly more compact form if we use the scaled 
coordinates 

p=— , C = -, a = -^. (45) 

Wo Zq 1 + C 

Then, the simplest wave function is 
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recalling eq. Q and the definition of /(C) in eq. (j32j) . In this manner the general, paraxial 
wave function can be written 

^n{r±, (p, z, t) = /'»+ip"^e--^'''L;7(2a)e*(^^-'^*±"''^-2"t^'^"' (47) 

It is noteworthy that although our solution began with the hypothesis of separation 
of variables in oblate spheroidal coordinates, we have found wave functions that contain 
the factors e"'^^ and L^{(t) that are nonseparable functions of r±_ and z in cylindrical 
coordinates. 

The wave functions found above are for a pure frequency u. In practice one is often inter- 
ested in pulses of characteristic width r whose frequency spectrum is centered on frequency 
u. In this case we can replace the factor e*^'^^"'^*) in the wave function by g{ip)e^'^, where the 
phase is (p = kz — ujt, and still satisfy the wave equation (jS} provided that the modulation 
factor g obeys [TT] 

q' 

- < 1. (48) 
9 

An important example of a pulse shape that satisfies eq. (jlH)) is 



giif) = sech^ , (49) 

UJT 

SO long as lot ^ 1, i.e., so long as the pulse is longer that a few periods of the carrier 
wave. Perhaps surprisingly, a Gaussian temporal profile is not consistent with condition 
Hence, a "Gaussian beam" can have a Gaussian transverse profile, but not a Gaussian 
longitudinal profile as well. 



2.2 Electric and Magnetic Fields of Gaussian Beams 

The scalar wave functions ()47j) can be used to generate vector electromagnetic fields that 
satisfy Maxwell's equations. For this, we use eqs. (fTSj) with a vector potential A whose 
Cartesian components are one or more of the functions ()47p. 

If we wish to express the electromagnetic fields in cylindrical coordinates, then we im- 
mediately obtain one family of fields from the vector potential 

A, = Ay = A^ = A^ = 0, A, = 0, ^, t). (50) 

The resulting magnetic field has no z component, so we may call these transverse magnetic 
(TM) waves. If index m = then A has no dependence, and the magnetic field has no 
radial component; the magnetic field lines are circles about the z axis. 

The lowest-order TM mode, corresponding to indices m = n = 0, has field components 

= Eopfe-fp"ge'^ + 0{el), 

= 0, 

= teoE,f\l-fp')e-fp'ge'^ + 0{el). 

B± = 0, 

Btj, = E±, 

= 0, (51) 
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as apparently first deduced in |221- This is a so-called axicon mode jTTj, in which the electric 
field is dominantly radial, which component necessarily vanishes along the beam axis. In 
the far zone the beam intensity is largest on a cone of half angle 6q and is very small on the 
axis itself; the beam appears to have a hole in the center. 

The radial and longitudinal electric field of the TMq mode are illustrated in Figs. El and 
EJ Photographs of Gaussian-Laguerre laser modes from [221 shown in Fig. |31 




Figure 2: The electric field Er{r±,0, z) of the TMq axicon Gaussian beam 
with diffraction angle 9o = 0.45, according to eq. (|KT|l . 



As is well known, corresponding to each TM wave solution to Maxwell's equations in free 
space, there is a TE (transverse electric) mode obtained by the duality transformation 

Ete = Btm, Bte = — Etm- (52) 

Since we are considering waves in free space where V • E = 0, the electric field could also be 
deduced from a vector potential, and the magnetic field from the electric field, according to 
the dual of eq. (fT^ . 

E = V X A, B = -|v X E. (53) 

k 

Then, the TE modes can be obtained by use of the vector potential in eq. (j^ . 
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ID 

Figure 3: The electric field Ez{r±,0, z) of the TMq axicon Gaussian beam 
with diffraction angle 9o = 0.45, according to eq. (jHTjl . 




Figure 4: Photographs of Gaussian-Laguerre laser beams with 2n + m = 2. 
From [23. 



The TM Gaussian-Laguerre modes emphasize radial polarization of the electric field, 
and the TE modes emphasize circular polarization. In many physical applications, linear 
polarization is more natural, for which the modes are well-described by Gaussian-Hermite 
wave functions [21 El IH El- Formal transformations between the Gaussian-Hermite wave 
functions and the Gaussian-Laguerre functions have been described in |24j . 
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2.3 Energy, Momentum and Angular Momentum in the Far Zone 

The electromagnetic field energy density, 

«=^^. (54) 

the field momentum density, 

E X B 

and the field angular momentum density, 

1 = r X p, (56) 

are the same for a TM Gaussian-Laguerre mode and the TE mode related to it by the duality 
transformation ()52j) . 

We consider the energy, momentum and angular momentum for TM waves in the far 
zone, where C ~ r/^o ^ 1, and r± ^ r6 <^ r in terms of spherical coordinates (r, 6, 0). Then 
the waves are nearly spherical, and so have a phase factor e*'^'" that implies the electric field 
is related to the magnetic field by 

E=-^VxB^Bxr, (57) 

rv 

SO that = B^. The time- averaged densities can therefore be written 

(^) = ^ , (58) 
on 

_ Re[{B X r) x B^] _ r - Re[{B ■ r)B-] ^ ^^^^ 
Stic Sire 

fle|(B.r)(.-xB')| 

^ ' 87TC ^ ' 

The TM waves are derived from the vector potential ()5()|1 whose only nonzero component 
is Az- Then, the magnetic field components in cylindrical coordinates are 

B± = TT- - , B^ - -- — — , i^^ - U, (bij 

r_L 0(p r± or± wq op 

where p = r±/wo. The radius vector r has cylindrical components (r^, 0, z), so 

B-r=^-^, and r x B* = -r^zB^ + 4)zBl + zr^B^. (62) 

Only the z component of the angular momentum can be nonzero for the beam as a whole, 
so we calculate 

(h) = --^ReiB^B;) ^ ORe lA,—^ = jr^Re lA,—^ 

Sircr ^ Sttcwq \ op J Sircwo kwqUq \ Op J 

- -^-Re[%A,^\, (63) 



^.-kujwq Oq \ dp 
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recalling that r^/r 9 and 9q = 2/ kwQ . The factors of that depend on r± and a ~ z^p^ /r"^ 



are 



A, oc p^e-f'''LZ{2a) ^ p^e-^^' ( 1 - ^^^'[yA , (64) 



where in the far zone, /(C) ~ —iz^jr. Thus, 

dAy^ I m Zn Anzn \ . Zn , 6 , , ^ 

— ^ ^ - + 22^ - ; pA, ^ 22^pA, ^ 2^-A„ (65) 

since in the far zone the factor p^e^^P^ ~ p™e~''^/^ implies that the wave functions are large 
only for p f» ^/m C r/2;o ^ 1- Inserting this in eq. we find 

To compare with the energy density, we need 

dA 





I 7-» 1 2 I 7-» 1 2 I T-» 1 2 I /I 1 2 

\B\ = \B_i\ + -Dfli =—5-/12 H 5- 

r1 Wn 



~ — — + 4^ 2 ~ ^ , 67 



9p 

since p = tttOq/XO ^ 9q/9 m the far zone. Using this in eq. (j58p . we have 

1 6^ 



and the ratio of angular momentum density to energy density of a Gaussian-Laguerre mode 
is 



(69) 



(n) LU 

where the ± sign corresponds to azimuthal dependence e^*"*"'^. 

In a quantum view, the mode contains photons per unit volume of energy hu each, so 
the classical result (j^Hl) implies that each of these photons carries orbital angular momentum 
±m^. Since the photons have intrinsic spin S* = 1, with 5"^ = ±1, we infer that the photons 
of a Gaussian-Laguerre mode carry total angular momentum component = ±m ± 1. 

The angular momentum of Gaussian-Laguerre modes has also been discussed in [22], in a 
slightly different approximation. The first macroscopic evidence for the angular momentum 
of light appears to have been given in fl^ . 

Using the above relations we can evaluate to momentum density (jSHI); which is propor- 
tional to the Poynting vector, and in the far zone we find 

Since kr± 3> 1 in the far zone, the energy flow is largely radial outward from the focal region. 
The small azimuthal component causes the lines of energy flow to become spirals, which lie 
on cones of constant polar angle 6 in the far zone 

We can, of course, also deduce the angular momentum density (l^) = [r x (p)]^ using 
eq. dZOD. 
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